Structure of end states for a Haldane Spin Chain 
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Inelastic neutron scattering was used to probe edge states in a quantum spin liquid. The ex- 
periment was performed on finite length antiferromagnetic spin-1 chains in Y2BaNii_xMgx05. At 
finite fields, there is a Zeeman resonance below the Haldane gap. The wave vector dependence 
of its intensity provides direct evidence for staggered magnetization at chain ends, which decays 
exponentially towards the bulk = 8(1) at T = O.I K). Continuum contributions to the chain end 
spectrum indicate inter-chain-segment interactions. We also observe a finite size blue shift of the 
Haldane gap. 

PACS numbers: 75.25.-l-z, 75.10.Jm, 75.40.Gb 



Highly-correlated electron systems can have quan- 
tum ground states that are coherent over macroscopic 
length scales. Among the best known examples is the 
two-dimensional electron gas in a magnetic field where 
Coulomb interactions lead to cooperative excitations 
with fractional charge as evidenced in the quantum Hall 
effect. Spin systems with strong quantum fluctuations 
feature similar cooperative phenomena with coherent 
spin excitations that differ fundamentally from classical 
spin waves. 

Impurities and boundaries can liberate novel compos- 
ite degrees of freedom and dramatically change the phys- 
ical properties of macroscopic quantum systems. Spin 
correlations in the vicinity of impurities are of particular 
interest because they reflect the nature of the underlying 
bulk quantum state, and because they may be important 
for understanding unconventional bulk properties. Many 
recent experiments have therefore been carried out to 
probe the effects of impurities on quantum fluids. For 
example, scanning tunnelling microscopy and NMR have 
examined the response of high-temperature superconduc- 
tors to Zn, Ni and Li substituted for Cu atoms ||, 
while neutron scattering has revealed that holes create 
AF droplets with a central tt phase shift in the one- 
dimensional Haldane spin liquid 

In this Letter we present direct information 
about vacancy-induced edge states in the one- 
dimensional (ID) S=l Heisenberg antiferromagnet (AF) 
Y2BaNii_xMgx05 obtained through a novel inelastic 
neutron scattering technique. The relative simplicity 
of a ID 5=1 antiferromagnet allows straightforward 
interpretation of the data, and therefore can serve as a 
step toward understanding impurity-induced phenomena 
in more complicated systems. The Hamiltonian for an 



AF spin chain is 



(1) 



where J > is the nearest- neighbor exchange constant. 
For S=l the ground state of Eq. |l| is an isolated singlet 
in which correlations fall off exponentially and the exci- 
tation spectrum has an energy gap [Q . The ground state 
is similar to that of the Affleck-Kennedy-Leib-Tasaki 
(AKLT) model where the S=l states are described 
as two S—^ states and the ground state is made up of 
local AF valence bonds between S—^ states of nearest- 
neighbor sites. The first excited state corresponds to 
breaking a valence bond singlet, hence the gap. 

Non-magnetic impurities in S=l chains should pro- 
duce an ensemble of finite length chain segments with 
a singlet-triplet gap that increases with decreasing chain 
length Q as well as localized S=^ degrees of freedom at 
chain ends The strongest evidence for S—^ degrees of 
freedom from non-magnetic doping was obtained through 
ESR § and NMR § experiments. It is important to 
note, though, that in real S=l chains, we are dealing 
with effective S=^ degrees of freedom, which always oc- 
cur in S—0 and S—l pairs, a fact graphically illustrated 
by magnetic specific heat data 1^. Here we present in- 
elastic neutron scattering experiments which probe the 
microscopic structure of the chain end states through the 
wave- vector dependence of the Zeeman resonance. 

Y2BaNi05 has a body-centered orthorhombic struc- 
ture, space group Immm with lattice parameters a = 
3.76A, b = 5.76A and c = 11.32A The Ni^+ 

chains run along the a- axis and the Ni^^ions are sep- 
arated by a, so Q—{h,k, I)— {0.5,0,0) is the AF point, 
which is also referred to as the q = 7r-point in theoretical 
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work. The samples were grown by the traveUing-solvent 
floating-zone method [Tl| . The experiments were per- 
formed at the NIST Center for Neutron Research. 

The excitation spectrum of AF S=l chains near AF 
wave-vectors is dominated by well-defined triplet excita- 
tions with an energy gap A = 0.41050(2)J Open 
circles in Fig. [l] show the spectrum of the nominally pure 
sample, Y2BaNi05, at the AF point q = tt, measured at 
(0.5, 1.5,0). A non-magnetic background, determined by 
matching constant-Q scans with the sample offset by 20°, 
was subtracted from the data. Crystal field anisotropics 
in Y2BaNi05 split the Haldane triplet into three sepa- 
rate modes at 7.5(1), 8.6(1) and 9.6(1) meV for polar- 
izations parallel to the three principal orthorhombic axes 
[p2|. The experiment at (0.5,1.5,0) is least sensitive to 
5^(Q, Lu) because neutron scattering probes the dynamic 
spin correlation functions of spin components perpendic- 
ular to the wave- vector transfer Q: 
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t^(|i^Ni(Q))'^(l-QaQ/3)^"''(Q,^). 

op 
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Here FNi(Q) is the magnetic form factor of Ni^+. 

Closed circles in Fig. |^ show the spectrum for the 
doped sample, Y2BaNio.96Mgo.04O5, measured at the 
same wave-vector and under similar resolution condi- 
tions. In agreement with work on powder samples of 
Y2BaNio.96Zno.04O5 [|l^, the spectrum is attenuated and 
harder than in the undoped sample, even though the 
threshold and peak energies remain essentially the same. 
This follows because there is now a distribution of chain 
lengths with the Haldane gap rising with decreasing chain 
length. The longest chains contribute unchanged spectra 
and threshold energies, while the shorter chain spectra 
are (blue) shifted to higher energies and produce the en- 
hanced - relative to the peak - high frequency tail. 

Indeed the observed neutron spectra follow quantita- 
tively from the random distribution of the non-magnetic 
defects in the sample. For defect density, x, the prob- 
ability that L + 2 sites form an isolated chain segment 
of length L is P{L,x) = x^il — x)'" . Numerical calcu- 
lations suggest that the energy gap in a chain of length 

L is A"(L) = ^ (A^)2 + 1,2 . sin2(7r(l - l/L)) §, where 
Vs = 2. 49(1) J is the spin-wave velocity [Q. The average 
dynamical structure factor is then given by S'(Q,cj) — 
P(L)S'i(Q, w), where S'l(Q,w) is the structure fac- 
tor for a chain segment of length L. For simplicity we 
chose Sl{Q,uj) to equal that for an infinite length chain 
suitably normalized so as to correspond to L spin and 
with gap A(i) I p^. ^(Q, to) was convoluted with the res- 
olution function |[16[ and adjusted to the data in a global 
fit. The fit for the nominally pure sample included a 
residual impurity concentration, x ~ 0.010(1), that was 
determined from the Curie tail in low-T susceptibility 
data p2] . For the doped sample the impurity concentra- 
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FIG. 1: Scattering spectra at (0.5, 1.5,0) for doped and un- 
doped Y2BaNi05, shown in absolute units per Ni^^. I{Q,lu) 
is the convolution of the cross-section with the resolution 
function. The solid lines are the scattering expected for the 
given distribution of chain lengths and convoluted with the 
resolution function. The nominally pure sample, Y2BaNi05, 
(m = 0.7 g) was studied using BT4 with 40'-40'-40'-80' col- 
limation, a final energy Ef = 14.7 meV and a PG filter in 
the scattered beam. The experiments on the doped sample, 
Y2BaNio.96Mgo.04O5, (m = 0.6 g) were performed using BT2 
with the coUimation set to 60'-40'-40'-80', the final energy 
Ef = 13.7 meV and a PG filter after the sample. Intensities 
were normalized using the integrated intensity of a transverse 
acoustic phonon close to the (200) reciprocal lattice point. 



tion was fixed at x = 0.04. 

The only adjustable parameter was a common overall 
scale factor, which in the single mode approximation is 
proportional to the ground state energy. The Haldane 
gaps A^ and the mode velocity Vg were fixed at the val- 
ues measured in the undoped compound [ p^ . This yields 
{Ti)/L ~ —1.45(0. 4) J, consistent with numerical results 
Q. The excellent agreement between experiment and 
model in Fig. |l| validates the assumptions made above. In 
particular the model accurately reproduces the increase 
of the average gap with impurity concentration x - and 
thus provides support for upward renormalization of the 
gap in accordance with A(L). 

We have thus far concentrated on the effects of impu- 
rities at energies near the Haldane gap A of the pure 
compound. To introduce what occurs at much lower 
energies, we recall that beyond the triplet excitations 
above A, the AKLT model predicts the existence of a 
pair of localized composite S=^ degrees of freedom to- 
wards the ends of finite length S=l chains. These are 
coupled via an indirect exchange interaction K{L) — 
(-l)^0.8064Jexp(-(L-l)/O 0, so that a doped 5=1 
chain system can be described as a set of dimers with 
a broad distribution of ferro- and antiferromagnetic in- 
terdimer couplings. The ground state of a chain seg- 
ment is a singlet when L is even and a triplet when L 
is odd. For chains that are much longer that the Hal- 
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dane length the gap between these levels is much smaller 
than the Haldane gap. A magnetic field, if, splits the 
triplet state and inelastic neutron scattering can then be 
used to probe transitions between quantum states of the 
nanometer-scale chain segments. 

k(r.l.u.) 
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FIG. 2: Spectra at Q = (0.5, fc, 0) for three different fields 
along the c axis in absolute units (left axis). The solid line is 
the model described in the text in absolute units (right axis), 
scaled so that the model matches the scattering at zero field. 
The top scale shows the median wave-vector transfer in the 
transverse direction which was varied to maintain kf paral- 
lel to the chain direction for optimal focusing. The experi- 
ment was performed with beam divergences 40'-80'-52'-300' 
for each of the 11 analyzer blade channels, which were set to 
reflect Ef — 3.7 meV to the center of the detector. A BeO fil- 
ter rejected neutrons with energies higher than 3.7 meV from 
the detection system. Normalization of the intensity relied on 
the incoherent elastic scattering from the sample. 

The consideration of the coupled chain end states sug- 
gests searching for excitations with much lower ener- 
gies than the Haldane gap. While such excitations can- 
not exist in pure chains, they should appear in diluted 
chains at energies of order the indirect exchange inter- 
actions, and will yield an inelastic tail extending into 
the gap from Tiuj = 0, where the chain end states would 
reside if the couplings between them were to vanish. 
Fig. ^ shows that as expected, this phenomenon oc- 
curs for Y2BaNio.96Mgo.04O5. To obtain the spectrum 
displayed, scattering at (0.4, /c,0) and (0.6, fc,0) was av- 
eraged and then subtracted as background. For finite 
fields the spin fluctuations, including those hidden by 
elastic incoherent scattering should be Zeeman-shifted. 
Our data also agree with this expectation - the neutron 



spectra at Q = (0.5, k, 0) (Fig. ^ and c) show a well- 
defined excitation, centered at the Zeeman energy, with 
clear low energy sidebands. 

The data in Fig. | can largely be explained using a 
model of uncoupled dimers consisting of the composite 
S=^ degrees of freedom at the chain ends. To include the 
spin anisotropics in the system we calculated the neutron 
spectra using the low-energy Hamiltonian introduced by 
Batista et al.[^: 

Hl = El + [JaL + Dl3L] \0){0\ 
+jLDS^,+jLE{Sl + Sl)-fiBgzH,S,, (3) 

where El, a^, and were computed by the density- 
matrix renormalization group technique and |0)(0| is 
the projection onto the singlet state. S is the spin oper- 
ator for the dimer. The effect of the anisotropy parame- 
ters D = -0.039 J and E = -0.0127J [|| is to split the 
triplet energy levels. 

The neutron scattering cross-section for a dimer is 

J^r(Q)^i(Q)EA.A',PA.(AL|exp(-iQ . Ri)5f |A'J 
(A'^l exp(zQ • Ri)S^\XL)S{Ex, - Ey^ + fiu) , (4) 

where denotes the dimer singlet and triplet states 
with energy E\^ and population factor . The in- 
dex i,j — 1,2 numbers the composite 5*—^ at a posi- 
tion Ri towards the end of the chain and a,(3 — x,y, z. 
We used g — = 2.165 as determined from ESR |l^ 
and (7ro)^ = 0.29 barn. -F2(Q) is the form factor of 
the end-chain composite spin and describes the decaying 
AF magnetization at the chain ends associated with the 
collective spin degree of freedom. We will give an ex- 
perimental view of this quantity below; for the moment, 
we simply assert that its square is well described by a 
normalized Lorentzian centered aX h = 0.5 with a width 
r given by = (^ + ■p-)/"'^. This should be a good 
approximation because the Fourier transform of an expo- 
nentially decaying correlation function for the staggered 
magnetization is a Lorentzian, and the finite length L of 
the chain acts as a finite size cutoff. 

The scattering cross-section for a given dimer distri- 
bution was obtained by adding contributions from all 
chain segment lengths -£§^ = Y.LPiL,x)-^{L) and 
convoluting with a Gaussian to account for the finite en- 
ergy resolution. Lacking detailed singlet and triplet wave 
functions for a chain segment, interference effects, such 
as the twin incommensurate peaks for the subgap bound 
states in hole-doped Y2_xCaxBaNi05 between chain 
ends were ignored. This should be a good approxima- 
tion for an ensemble of long chains. The solid line in 
Fig. H shows the absolute scattering intensity expected 
from this model. The weakly oscillatory nature of the cal- 
culated spectrum follows from the discrete nature of the 
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distribution of indirect exchange couplings K{L). At zero 
field there is excellent agreement between experiment and 
theory with an 28% adjustment of the pre-factor (within 
the 30% error bar from normalization) being the only ad- 
justable parameter. At finite field the model reproduces 
the position of the well-defined peak {g^ — 2.0(0.05) gives 
a sfightly better fit than = 2.165 from ESR but 
it overestimates its intensity, which seems to be trans- 
ferred to a continuum that cannot be accounted for by 
the model. The broad scattering below the well-defined 
excitation indicates that changes in the chain end struc- 
ture with field or interchain coupling play an important 
role. If interchain interactions indeed couple the chain 
segments then the doped sample may eventually freeze 
into a low-temperature AF or spin glass state. 

Neutron scattering has the tremendous advantage of 
providing spatial as well as spectral information about 
spin fluctuations. We have thus examined the form fac- 
tor for the edge states at the Zeeman energy by scanning 
wave vector h with fixed huj. Fig. |^ shows the result. 
Rather than a broad maximum centered at h — corre- 
sponding to the form factor for a localized atomic spin, 
we find a relatively sharp peak centered at h = 0.5. This 
is direct evidence that chain end spins carry AF spin po- 
larization back into the bulk of the chain segments. By 
fitting the model cross-section (Eq. ||) convoluted with 
the resolution function to the data we find the exponen- 
tial decay length ^ for the chain-end staggered spin den- 
sity (inset of Fig. |^). The decay length remains constant 
at ^ — 8(1) between H — 4 and 11 T - in perfect agree- 
ment with the equal-time correlation length ^ = 8(1) 
for undoped Y2BaNi05 in zero field at low temperatures 
T/J 1 |l^, but slightly longer than that deduced from 
an NMR study of the staggered magnetization 

[§• A fit 

to a square-root Lorentzian, which would be appropri- 
ate for unpinned (by impurities) bulk spin correlations 
from the Haldane continuum, yields a longer correlation 
length ^ — 13(1). Note that because we have fixed hw 
at precisely the Zeeman energy, we have maximized the 
contribution of those chains with the smallest coupling 
and hence largest lengths to the scattering, thus reveal- 
ing F£^^(Q), in the notation of Eq. ^. 

In summary, our inelastic neutron scattering experi- 
ment on Mg-doped Y2BaNi05 shows an upward renor- 
malization of the Haldane gap with doping, which is con- 
sistent with numerical predictions for a distribution of fi- 
nite chain segments. In addition, we find well-defined ex- 
citations below the Haldane gap whose field dependence 
demonstrates that they originate from 5=^ chain-end 
degrees of freedom. Their wave- vector dependence shows 
chain end states carry an AF "Friedel" oscillation that 
extends into the chain segment with an exponential decay 
length — 8(1) in a field H = IIT. The experiment also 
demonstrates a new technique, neutron ESR, by which 
the structure of composite spin degrees of freedom liber- 
ated from a macroscopic singlet can be probed. Efforts 
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FIG. 3: Scattering originating from S—^ chain end states at 
the resonance energy huj = 1.30 meV for H — 11 T (applied 
along the c axis) in Y2BaNio.96Mgo.a4O5 as a function of h. 
The solid line is a fit explained in the text and the dashed line 
is the background. Inset: Correlation length ^ as a function 
of field. For this experiment Ef = 5 meV and there was a 
cooled Be filter in the scattered beam. 



are now underway to apply this approach to understand 
impurity structure and dynamics in more complex sys- 
tems such as high-temperature superconductors pc| ] . 
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